function [population] = InitializePopulation(Dimension,Particle_Number,VRmin,VRmax)

ps=Particle_Number;
D=Dimension;
rand('state',sum(100*clock));
if length(VRmin)==1
    VRmin=repmat(VRmin,1,D);
    VRmax=repmat(VRmax,1,D);
end
mv=0.5*(VRmax-VRmin);
VRmin=repmat(VRmin,ps,1);
VRmax=repmat(VRmax,ps,1);
Vmin=repmat(-mv,ps,1);
Vmax=-Vmin;
population=VRmin+(VRmax-VRmin).*rand(ps,D);

end